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Abstract 

We investigate in this paper the antideuteron flux produced in high energy colHsions 
of cosmic rays with the interstellar matter. We employ the Monte Carlo generator 
DPMJET-III together with the coalescence model to simulate, in an event-by-event 
basis, the antideuteron production in cosmic ray collisions in our Galaxy. Then, we use 
a diffusion model to calculate the expected flux. We find a secondary antideuteron 
flux at the Earth which is approximately a factor of two smaller than in previous 
calculations. 

1 Introduction 

The existence of cosmic antiprotons produced in high energy collisions of cosmic rays with 
the interstellar gas was suggested by Fradkin P shortly after the discovery of the antiproton 
at the Bevatron [2] and positively observed as a cosmic ray component in the late seventies [3l 
H]. The measurement of the spectrum was refined in subsequent experiments O El [71 
[HI [9], culminating with the exquisite recent measurement of the antiproton flux and the 
antiproton-to-proton fraction by the satellite borne experiment PAMELA [TT] . 

The discovery of antinucleons opened the possibility of observing antinuclei in Nature. 
Indeed, antideuterons were first observed in 1965 in proton-beryllium collisions at the CERN 
proton-synchrotron [12]. It is then natural to expect the observation of antideuterons in 
cosmic rays, although so far all searches have been fruitless. The best present limit on 
the cosmic antideuteron flux was set by BESS in the range of kinetic energy per nucleon 
0.17 < T < 1.15GeV/n, $j < 1.9 x 10"^ m-2s-isr"i(GeV/n)"^ [I3]. Interestingly, in the 
near future, the sensitivity of experiments to the cosmic antideuteron flux will increase sig- 
nificantly, by more than two orders of magnitude. The Alpha Magnetic Spectrometer Exper- 
iment (AMS-02) on board of the international space station is currently searching for cosmic 
antideuterons in two energy windows, 0.2 < T < 0.8 GeV/n and 2.2 < T < 4.4 GeV/n, 
with an expected flux sensitivity after five years $j = 1 x 10^^m~^s^^sr^^(GeV/n)~^ in 
both energy windows [11]. Furthermore, the balloon borne General Antiparticle Spec- 
trometer (GAPS) will undertake, starting in 2014, a series of flights at high altitude in 
the Antarctica searching for cosmic antideuterons. In the first phase, a long duration 
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balloon (LDB) flight will search for antideuterons in the range of kinetic energy per nu- 
cleon 0.1 < T < 0.25 GeV/n, with a sensitivity $j = 1.2 x 10-^m-^s-^ST-^{GeV /n)~^ , 
while in the second, the ultra long duration balloon (ULDB) flight will search in the range 
0.1 < T < 0.25 GeV/n, with a sensitivity $j = 3.5 x 10-^m-2s-isr-i(GeV/n)"^ [H]. 

The expected fluxes of cosmic antideuterons produced in collisions of high energy cos- 
mic rays with the interstellar gas have been previously investigated in flbl 116] . A precise 
calculation of the secondary antideuteron flux is important not only to assess the prospects 
to observe antideuterons at AMS-02 or GAPS, but also to determine the background fluxes 
to search for primary antideuterons from exotic sources, such as dark matter annihila- 
tions dZl HSl 1121 EOl El ES], dark matter decays [211201 ED ES] or black hole evapora- 
tion [25]. Unfortunately, the calculation of the secondary fluxes is hindered by the scarce 
experimental information on the differential cross section for antideuteron production in high 
energy collisions, hence many approximations and ad hoc assumptions had to be adopted. 

In this work we aim to circumvent some of these assumptions by employing a Monte Carlo 
generator together with the coalescence model to simulate in an event-by-event analysis 
the antideuteron production in high energy proton-proton and antiproton-proton collisions. 
Our method to calculate the source term of secondary antideuterons is described in Section 
[2j We then calculate in Section [3] the antideuteron flux at the Earth using a a diffusion 
model, including convection, energy loss due to ionization. Coulomb collisions in a ionized 
plasma and non-annhilating rescatterings on the interstellar medium, as well as diffusive 
reacceleration caused by interactions with random magnetohydrodynamic waves. Finally, 
we compare our predicted fluxes to the results of previous literature and to the expected 
sensitivity of the experiments AMS-02 and GAPS, and we discuss the impact of the various 
sources of uncertainty in our flnal result. Our conclusions are presented in Section |4j 

2 The source term of secondary antideuterons 

The source term of secondary antideuterons at the position r with respect to the Milky 
Way center, Q^'^'^ [T^, r) , deflned as the number of secondary antideuterons produced per 
unit volume, time and kinetic energy per nucleon Tj, reads [T3] 



In this expression, i runs over all relevant incident cosmic rays species, with flux (Tj, r) and 
which are known at the Earth with fairly high accuracy. Besides, j represents the different 
components of the interstellar medium (ISM) which we assume uniformly distributed over 
the galactic disk, which extends radially 20 kpc from the center and has a half-thickness h of 
100 pc. Lastly, ^""''jj^^^d^ denotes the differential cross section to produce an antideuteron 
with kinetic energy per nucleon in a collision of a cosmic ray particle of type i (with 
kinetic energy per nucleon Tj) with an ISM component of type j. 

To describe the antideuteron production we will employ the coalescence model [201 EZl 
[28| l29l [20] . which postulates that the probability of formation of one antideuteron out of 
an antiproton-antineutron pair with given four-momenta kp and can be approximated 
as a narrow step function 9 (A^ -|- pi), where A'^ = kp — k'^. In this model, the coalescence 
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momentum pq is the maximal relative momentum of the two antinucleons that still allows 

the formation of an antideuteron. One can show |20j that for Ifc^l ^ po, being 

this ansatz leads to the following differential antideuteron yield in momentum space: 

d^N.^ 143 d'Npd'Nn ( h h\ 

^'IPk,^^'^ = 8 ■ ■ ^^^-l^%d^ U' ^ j ■ ^ ^ 

In absence of a microscopic understanding of the coalescence mechanism, the coalescence 
momentum po should be determined from experiments. It is important to emphasize that 
the coalescence momentum is not a universal parameter, but it shows a dependence on the 
underlying process and also on the center of mass energy of the corresponding reaction [23] . 
Therefore, the coalescence momentum employed for the antideuteron production in spal- 
lations of cosmic rays on the interstellar medium should be extracted from a laboratory 
experiment of antideuteron production in nucleus-nucleus collisions at the corresponding 
energy. 

It is common in the literature to assume the statistical independence of the antiproton 
and antineutron production in Eq.([2]). In this case the antideuteron yield is proportional to 
the antiproton yield times the antineutron yield. This assumption is, however, not justified 
for the production of low energy antideuterons in cosmic ray collisions, which is the energy 
range of interest for experimental searches, and which takes place very close to the threshold. 
More specifically, a pp collision at -y/s = 10 GeV producing one antiproton in the final state 
must contain at least three protons, reducing the available energy for other particles in the 
final state to ~ 6 GeV or less. The production of an antideuteron requires an additional 
antineutron and neutron, which is further supressed by the limited phase space available 
for these particles. This leads to a correlation of p and n production in collisions close to 
the antideuteron production threshold A/S(thres) ~ 6mp: having already one antinucleon in 
the final state, the probability of creating another one is reduced. Therefore the factorized 
coalescence model, which assumes uncorrelated production of the antinuclei, overpredicts 
the antideuteron yield. 

This drawback of the factorized coalescence model was circumvented in [TB] by intro- 
ducing a phase space suppression factor in Eq.([2]). Defining the relativistic phase space of 
n particles with masses rrii as 

$ (y^; mi, m2, . . . , m„) = 1 n ) (T^P^ ^ ' v^) ' 

and introducing a phase space suppression factor _R„ (x) for n nucleons with masses nip 
through 

^"^"^ = $(x;0,0,...,0) ' 

the "modified factorized coalescence model" states that the differential antideuteron yield 
in momentum space reads: 
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Figure 1: Antideuterons from T decay. Both calculations are using po = 133 MeV; see the 
text for details. 



where ^/s is the total available center-of-mass energy and m;i{Eg) is the mass (energy) of 
the antideuteron. The index n of the phase space suppression factor is the minimal number 
of baryons produced together with the antideuteron and reads n = 4 for a, pp collision and 
n = 2 for pp collision. The modified factorized coalescence model then takes into account 
the correlation between the antiproton and the antineutron in the pp collision due to the 
limitied phase space. However, it does not take into account possible correlations in the 
underlying hard process. The importance of these correlations is manifest in FigjT| where 
we compare the spectrum of antideuterons from T (IS) decays predicted by an event-by- 
event analysis implemented in PYTHIA 8 [30] for the value po = 133 MeV (taken from [23]) 
with the spectrum predicted by the modified factorized coalescence model, calculated using 
Eq. ^ with an antiproton spectrum simulated also with PYTHIA 8 and for the same value 
of the coalescence momentum po = 133 MeV. We also show for comparison the spectrum 
measured by CLEO [31]. As apparent from the plot, taking into account the correlations in 
the hard process reduces the antideuteron yield by approximately a factor of two. 

In view of this limitation of the modified factorized coalescence model we will employ 
in this paper a full event-by-event simulation using the event generator DPM JET-HI ^32], 
an implementation of the two-component Dual Parton Model that is in particular designed 
for low-energetic hadron-hadron collisions, such as those relevant for experiments searching 
for cosmic antideuterons. Despite the fairly accurate results of DPMJET-III in reproducing 
the average multiplicity of protons, pions and kaons in pp collisions [33], DPMJET-III 
overproduces antiprotons in the most relevant region of center of mass energy. This is 
illustrated in Fig. [2| left plot, which shows the antiproton production differential cross 
section in pp collisions at an incident kinetic energy Tp = 18.3 GeV predicted by DPMJET- 
III compared to experimental data from the S61 experiment at the CERN PS [3^. Evidently, 
DPMJET-III overestimates the p production cross section by a more or less constant factor 
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Figure 2: Left plot: Antiproton production cross section in pp collisions predicted by 
DPM JET-III compared to the experimental data from [51] , where 9 is the production angle 
of the antiproton, while the antiproton momentum kp is fixed to 4.5 GeV/c. We also present 
the corresponding predictions scaled by the function 5'(Tp), as defined in the text. Right 
plot: Scaling function S (Tp) as a function of the proton kinetic energy. 



of ~ 2. Therefore we use DPMJET-III in a slightly modified form. First, we extract the total 
antiproton multiplicity per pp collision as a function of the kinetic energy of the impinging 
proton from DPMJET-III, n?^^'^^"'" (Tp), and compare it to an interpolating function of 
experimental values Tip"^^ (Tp) [35] . Then we define a scaling factor via 
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This function parametrizes the overproduction of antiprotons in DPMJET-III and is 
shown in Fig. [2| right plot. The behaviour for Tp > 1 TeV, where no experimental informa- 
tion on the antiproton multiplicity is available, is irrelevant due to the small contribution 
of these processes to the antideuteron source spectrum; for definiteness, we choose S (Tp) to 
be constant for Tp > 1 TeV. The antideuteron production cross section in pp collisions is 
finally obtained via 
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where the cross section on the right hand side is evaluated implementing the coalescence 
model, Eq. (|2]), in DPMJET-III with a coalescence momentum pq = 152 MeV, deter- 
mined in [2JJ from the measurement of the antideuteron production data in pp collisions 
at = 53 GeV at the CERN ISR [Ml EZj- The ansatz Eq. uses the well-motivated 
assumption that the antineutron and antiproton spectra are equal, and furthermore that the 
antideuteron yield is proportional to the product of the antiproton and antineutron yield. 

In this way, we calculate for 53 logarithmically binned values of the incident kinetic 
energy Tp in the range between 25 GeV and 10 TeV, which suffices our purposes. These cross 
sections can finally be inserted in Eq. (jT| in order to obtain the pp contribution Q^^ (Tj) to 
the secondary antideuteron source spectrum. For this, we furthermore use nn = 1.0 cm~'^ 
and take the interstellar proton flux $p measured by the AMS-01 measurement [38] . 
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Figure 3: Contributions of different ranges of incident kinetic energies Tp to the antideuteron 
source spectrum of tfie pp cfiannel (left plot) and the pp channel (right plot). For the pp 
(pp) channel, (1) refers to the range 25-42 GeV (7-10.6 GeV), (2) to 42-94 GeV (10.6-27 
GeV) , (3) to 94 - 473 GeV (27 - 67 GeV) and (4) to 473 GeV - 10 TeV (67 - 1780 GeV). 
The sum of all contributions is shown in black. 



In Fig. |3| left plot, we show the secondary antideuteron source spectrum Q^^^ (Tj) origi- 
nating from collisions of cosmic ray protons on the interstellar hydrogen (black curve), to- 
gether with the different contributions from various ranges of incident proton kinetic energies 
Tp. It is apparent from the plot that the production of antideuterons with Tj < 1 GeV/n is 
highly suppressed, due to the large kinetic energy required for an incident proton to produce 
an antideuteron, T^^^^^^ = 16mp. As a result, all produced particles are boosted relative to 
the center-of-mass frame in the direction of the incident proton, thus being highly unlikely to 
produce an antideuteron that is almost at rest in the galactic frame. More specifically, one 
can show that the smallest possible kinetic energy per nucleon of a produced antideuteron is 
T™™ = rrip/A = 0.23 GeV/n, which requires a proton kinetic energy Tp > 1 TeV. However, 
due to the fast decrease of the primary proton flux with the energy, oc E""^'^, the num- 
ber of low energy antideuterons produced by interactions of high energy protons with the 
interstellar gas is accordingly very small. 

High energy interactions of cosmic ray antiprotons on the interstellar hydrogen and 
helium can lead, as first noticed in [15], to a non-negligible contribution to the secondary 
antideuteron source spectrum. In this process, there are only two additional baryons in 
the final state, hence the threshold for antideuteron production in a pp reaction is only 
rpthves ^ Q^^^ ^Yiile it is Tp*'^'"^^ = 16mp for the pp process. Therefore, in the energy range 
relevant for searches for cosmic antideuterons the production cross section in the pp channel 
can be larger than in the pp channel, partially compensating the smaller incident cosmic 
ray flux of antiprotons in comparison to protons. 

The calculation of the antideuteron yield in the process pp ^ dX within the modified 
coalescence model requires the knowledge of the differential cross sections for pp ^ pX and 
pp ^ nX. Unfortunately the differential cross sections for these two processes are largely 
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unknown, hence previous calculations [151 [IB] assumed 

d"^(j d^(j 

iPP ^ - Ep-^_ (pp ^ pX) , (7) 

where the differential cross section for pp pX is well known experimentally (see e.g. 
Ho]). However, the approximation ([t]), when evaluated at the relevant energy range ^/s 
10 GeV can be questioned from purely kinematical grounds using baryon number conserva- 
tion: in the reaction on the left hand side of ([T]), at least one additional baryon has to be 
produced, while the reaction on the right hand side has to have at least three additional 
baryons in the final state. Therefore one can expect that the latter reaction is suppressed 
in comparison to the former and hence this approximation used for the evaluation of the 
antideuteron production is questionable. 

In order to circumvent the limitation of this approach we then simulate the reaction 
pp — 7- dX using also DPMJET-III on an event-by-event basis. First, we test the validity 
of DPMJET-III in pp collisions by comparing the predictions of the (anti)proton yield in 
the reaction pp pX at Tp = 31.1 GeV with the measurements at the Mirabelle bubble 
chamber at IHEP [H]. As can be seen in Fig. |4]the agreement is fairly good (within ~ 20%) 
without any rescaling of the (anti)proton yields. Therefore we use DPMJET-III without any 
modifications for performing the event-by-event analysis of antideuteron production in pp 
collisions. In the absence of any information on the coalescence momentum in this process, 
we use the same value po = 152 MeV as for the pp collisions. Following this procedure, we 
calculate for 49 logarithmically binned values of Tp between 7 and 1780 GeV. Lastly, 
we calculate the secondary source spectrum of antideuterons originated in pp collisions, 
Q'^ppiT^), using Eq. ([T]) and adopting the values from |12] for the interstellar antiproton flux 
$p. The result is shown in Fig. [3| right plot. The antideuteron spectrum extends to slightly 
lower values of Tj compared to the pp channel, due to the fact that the minimal kinetic 
energy of the incident particle is only Tp^^'^^ = Grrip and therefore the produced particles 
are less boosted on average. Note that due to the much smaller interstellar antiproton flux 
{p/p ~ 10~^), we have only considered incident antiproton energies up to only ~ 2TeV. 

Lastly, we also include in our calculation of the source term of secondary antideuterons 
the contributions from the channels p He, He p and p He. We obtain the spectrum of 
antideuterons produced in a reaction of two nuclei A and B from the corresponding an- 
tideuteron spectrum in pp collisions by scaling it with an appropriate nuclear enhance- 
ment factor which, following [53], we assume independent of the energy of the final 
antideuteron: 

^ (Ta, T,) = (Ta) ■^ATv = Ta, T,) , (9) 

T4 being the kinetic energy per nucleon of the nucleus A. We estimate e^g for the relevant 
processes using the "wounded nucleon" model |13]. With this, the antideuteron production 
cross sections in the reactions p He, He p and p He can be obtained using the corresponding 
nuclear enhancement factor together with the results of our event-by-event analysis of the pp 
01 pp reaction, respectively. The contributions to the secondary source spectrum are then 
calculated via Eq. ([T]), using nne = 0.07 cm~'^ and the helium flux measured by AMS-01 [38] . 
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Figure 4: Comparison of antiproton production cross section in pp collisions in DPMJET-III 
with experimental data from ^41j . 



The total secondary antideuteron source spectrum Q^'^'^ (Tj), i.e. the sum of the five 
production channels pp, pHe, Hep, pp and pHe is shown in Fig. [5] together with the in- 
dividual contributions. Evidently, for > IGeV/n the most important contribution is 
the antideuteron production in pp collisions, while for Tj < IGeV/n the antiproton in- 
duced channels dominate. As elaborated above, this can be explained through the different 
kinematics of the pp and pp reaction, concretely through the different thresholds for an- 
tideuteron production in these processes. The sum of the production channels involving 
helium either as a cosmic ray particle or as a component of the interstellar matter con- 
tributes on average ~ 30% to Q^^"^ (Tj). 



3 Propagation 



Antideuterons produced in collisions of cosmic rays in the interstellar medium propagate 
in a complicated way before reaching the Earth. The number density of antideuterons in 
interstellar space, / {E,r,t), satisfies the following propagation equation [HI HJl 

df{E,f,t) ^ 



dt 



K{E, f) V/ -VJ]+ 2h 6{z) g^DEE{E, r)—f - 2h 6{z)T^^J 



d_ 
dE 



2h5(z) ^ ^ 

Coul 

g^^^ (E, f, t) + Q'^^ {E, f, t) 



dE 
~dt 

tcr 



dE 
~dt 



Ion. 



/ 



E' 



m-, 



3E 



(10) 



where E is the total energy of the antideuteron. We also require the solution f{E,f,t) 
to vanish at the boundary of the diffusion zone, which is approximated by a cylinder with 
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Figure 5: Contributions to the secondary antideuteron source spectrum Q**^^ (Tj) 



half-height L = 1 — 15kpc and radius R = 20kpc. The antideuteron flux $ can then be 
derived from the number density using 

^(E,r,t) = ^f{E,r,t) , (11) 

with v{E) being the antideuteron velocity. 

We assume in the following that the number density of antideuterons is stationary, 
^ = 0, since the typical time scale of antideuteron propagation is small compared to the 
time scale on which galactic propagation conditions change |16] . The flrst term in the right- 
hand side of the propagation equation accounts for the diffusion of charged particles in the 
Milky Way as a result of scattering on inhomogeneities of the galactic magnetic flelds, leading 
to a random walk of the particles. In our calculation we assume that the corresponding 
diffusion coefficient is constant in the whole diffusion zon^ and is parametrized as [IHl HI] 
K{E) = Kof37l^, where (3 = v/c and 71 is the rigidity, defined as the momentum per unit 
charge. The second term describes the convective transport induced by a galactic wind 
of charged particles which we model as = VcSgn{z) Cz The third term accounts 

for the diffusive reacceleration in the galactic disk caused by interactions with random 
magnetohydrodynamic waves, commonly taken into account by diffusion in energy space. 
Following imilQ], we use a spatially constant energy diffusion coefficient Dee{E) which is 
related to the spatial diffusion coefficient K{E) by 

OME)-lviffy (12) 

being Va the Alfven velocity. The fourth term describes particle losses due to annihilations 
of antideuterons with protons or helium nuclei in the interstellar gas in the galactic disk 

^For a discussion of the effects of a position dependent diffusion coefficient see, e.g. [T7] . 
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with a rate Fann = (^h + 4^/'^nHe) cr'^ v, being tt-h = 1 cm~^ and riRe = 0.07 cm~'^. Here it is 
assumed that the annihilation cross section between an antideuteron and a hehum nucleus is 
related to the annihilation cross section between an antideuteron and a proton by a simple 
geometrical factor 4^/^. The annihilation cross section of an antideuteron with a proton 
is obtained substracting the non-annihilating rescattering (NAR) cross section from the 
inelastic cross section: ctt"'^ = o"?"^^ — a^^^. The inelastic dp cross section can be estimated 

dp dp dp ^ 

through cr™^' (Tj) ~ 2 0"™^^ {Tp = Tj) [16] , using the inelastic pp cross section from [51] : 



4^ (Tp) = 24.7 (1 + 0.584 Tr°-ii5 + 0.856 Tr^'^^^) . (13) 



Besides, for the non-annihilating rescattering cross section, dj^^ (T^): we use the estimation 
presented in [TBj based on experimental data on the total cross section for the charge sym- 
metric reaction p-f-d — )• p/n + d+{n7r). The fifth term describes continuous energy losses in 
the galactic disk, originating from Coulomb collisions in an ionized plasma and from ioniza- 
tion losses, where a part of the energy of the antideuteron is transferred to an electron of an 
interstellar matter atom or molecule. For the energy loss rates and (^)jq^ we refer 

to the formulas presented in [S2]- The sixth term describes adiabatic energy gains or losses 
which result from non-uniform convection velocities [53], which, adopting Vc = V^sgn(2;) e^, 
is also proportional to 6 (z) and hence confined to the Galactic disk. The last two terms 
account, respectively, for the antideuteron production in high energy interactions of cosmic 
rays in the galactic disk, being the source term of secondary antideuterons described in 
Section [2] , and for the antideuteron energy loss through non-annihilating rescattering on 
the interstellar matter in the galactic disk, i.e. reactions of the type d + p — >■ li-l-Xor 
(J + He — )■ d + X. This process, dubbed tertiary source term, does not change the total 
number of antideuterons but leads to a redistribution of the antideuteron number density 
towards lower energies. Formally, the tertiary source term reads |15j : 



da^NAR {d{TC) + p d{T,) + X) 



- a]^"" {d (Tj) + p -> J + X) (T,-, ^ 



(14) 



where $ j (Tj) is the antideuteron flux at kinetic energy per nucleon Tj and cr^ (x^^ jg 



the non-annihilating rescattering cross section, introduced after Eq. (13). For the NAR 
differential cross section there is unfortunately no experimental information, hence we cal- 
culate it following the approach pursued in [TSl dS] which is based on the measurement 
by Anderson et al. of the inclusive p + p p + X cross section at 10, 20 and 30 GeV/c 
incident proton momentum [M]. The tertiary contribution coming from the NAR process 
(i + He — )■ (i + X is also included by assuming that all relevant cross sections can be obtained 
from the corresponding process involving a proton by multiplication by a geometrical factor 
42/3 Notice that the tertiary source term depends itself on the antideuteron flux, hence 
the propagation equation is an integro-differential equation which we solve following the 
method presented in [19] . 

In the propagation equation there are five undetermined parameters: the half-height of 
the diffusion zone L, the normalization of the diffusion coefficient /Tq) the spectral index 5 of 
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Model 


5 


Ko (kpcVMyr) 


L (kpc) 


K (km/s) 


Va (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


22.4 


MED 


0.70 


0.0112 


4 


12 


52.9 


MAX 


0.46 


0.0765 


15 


5 


117.6 



Table 1: Definition of MIN, MED and MAX propagation parameters (see text). The values 
are taken from |55j . 



the rigidity in the diffusion coefficient, the convection velocity Vc and the Alfven velocity Va 
which have to be determined from observations. We list in Table [l] three sets of parameters 
derived in [55] compatible with the measurements of the Boron-to-Carbon ration and which 
yield the the minimal (MIN), medium (MED) and maximal (MAX) primary antiproton 
flux. The dependence of the final flux on the chosen set of parameters can therefore be 
considered as the uncertainty band connected to the imprecise knowledge of the propagation 
parameters. 

Lastly, in order to calculate the flux of a charged cosmic ray species at Earth, one has to 
take into account the effect of the solar wind. We adopt in this paper the effective theory 
of a spherically symmetric and charge-independent force field presented in [561 EZl EH] • In 
that simplified model, the net effect of the solar wind is an electric potential generated by 
the Sun leading to an energy loss of each charged particle. Then the flux at the top of the 
Earth's atmosphere •l)'^^'^ is related to the interstellar flux through 



$TOA (Ttoa) - ^$is (Tis) - 1^ 2mATis + (Aris)2 ' ^^''^ ^^^^ 



for a cosmic ray nucleus with mass m and mass number A. As usual, T denotes the kinetic 
energy per nucleon, i.e. T = ^j^- The top of atmosphere and interstellar kinetic energy 
per nucleon are related via Tjs = Ttoa + a'''^^^' "^^^^^ Z is the electric charge of the nucleus 
in units of e and 0f is the Fisk potential. The latter quantity is the effective parameter 
controlling solar modulation and varies between 500 MeV and 1.3 GeV over the eleven year 
solar cycle. In our calculations we will adopt the value 0f = 500 MeV (corresponding to 
solar minimum) and we will discuss later the impact on our results of using a different value 
of the Fisk potential. 

We show in Fig. [6] our final result for the antideuteron secondary flux at the top of the 
Earth's atmosphere (TOA), using a Fisk potential 0f = 500 MeV and the MED propaga- 
tion model. Note that the antideuteron background flux extends to arbitrary small kinetic 
energies, in contrast to the secondary source spectrum calculated in Section |2] and shown 
Fig. [sj which presents a hard cut-off at Tj ~ 0.3 GeV/n. This is due to the fact that for 
very low kinetic energies the background flux is dominated by the adiabatic energy loss 
and the tertiary effect, both significantly redistributing antideuterons from intermediate 
(Tj ~ 5 GeV/n) to low (Tj ~ 0.1 — 1 GeV/n) kinetic energies. We also show, for com- 
parison, the results for the TOA antideuteron secondary flux obtained in the most recent 
calculations [151 US] , which use the modified factorized coalescence model instead of a full 
event-by-event simulation, which results in different antideuteron production cross sections 
in the high-energy reactions. Lastly, we also show in the figure the prospected sensitivities 
of AMS-02 and GAPS ULDB. Both experiments are clearly not sensitive enough to de- 
tect antideuterons from collisions of cosmic rays in the interstellar matter. Translated into 
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Figure 6: Antideuteron background flux at the top of atmospfiere (TOA) for the MED 
propagation model and a fisk potential ipp = 500 MeV. The solid line corresponds to the 
result of this work, while the dotted and dashed lines are taken from Duperray et al. [15] and 
from Donato et al. [H] , respectively. The shaded regions indicate the prospected sensitivities 
of AMS-02 and GAPS ULDB. 



number of events, one expects 0.13 events at AMS-02 and 0.024 events at GAPS ULDB. 
Therefore, the observation of antideuteron events at AMS-02 or GAPS can be interpreted 
as a signal for an exotic source of antideuterons, although, as argued in [23], the stringent 
limits on the primary antiproton flux from the PAMELA measurements of the p/p ratio 
make this possibility also unlikely. The searches for cosmic antideuterons away for the solar 
minimum are more pesimistic, due to the larger solar modulation effects. This is illustrated 
in Fig. [Tj, where we show the top of the atmosphere antideuteron flux for values of the Fisk 
potential (pp = 500 MeV and (pp = 1.3 GeV, as well as the interstellar flux. 

The calculation of the cosmic antideuteron flux at the Earth suffers from a number of 
uncertainties. The largest uncertainty originates from the unknown value of the coalescence 
momentum po in the interactions of cosmic ray nuclei with the interstellar gas. For our 
calculation we use po = 152 MeV, which is the value we extract from pp collisions at 
a/s = 53 GeV using the coalescence model implemented in the event-by-event generator. 
Nevertheless, as discussed in Section [2] (see also [23j) the coalescence momentum displays 
a dependence on the underlying process and on the center of mass energy involved. We 
then investigate in Fig. [8| upper left plot, the impact of varying the value of the coalescence 
momentum on the cosmic antideuteron flux, showing the result for p^ = 133 MeV, as 
extracted from the T decay, and for pq = 192 MeV, as extracted from the Z decay [23] . 
We would like to stress here that this source of uncertainty arises due to the very scarce 
experimental information on antideuteron production in pp collision, which led us to use a 
Monte Carlo generator to simulate the production. This source of uncertainty could then 
be greatly reduced with a dedicated experiment measuring the differential antideuteron 
production cross section in pp collisions at y/s ~ 10 GeV (or Tp ~ 50 GeV in a fixed target 



12 




Figure 7: Antideuteron background flux in interstellar space (IS) and at the top of the 
Earth's atmosphere (TOA), for two different values of the Fisk potential (pp. All curves 
assume MED propagation parameters. 



experiment), which is the most relevant process and center of mass energy for the calculation 
of the cosmic antideuteron flux, cf. Fig. |3| 

The second largest uncertainty of our computation of the cosmic antideuteron flux stems 
from the degeneracies in the propagation parameters. We show in Fig. [8| upper right plot, 
the resulting flux for the three different sets of propagation parameters, MIN, MED and 
MAX, in Table [Tj As apparent from the plot, the MIN and MED propagation parameters 
lead to nearly the same result, while the MAX parameters lead to a cosmic flux which is 
lower by ~ 40 % in the low-energy region. The reason for this dependence is the following: 
for kinetic energies above ~ 5 GeV/n, where energy loss effects are unimportant, the impact 
of the choice of propagation parameters is rather mild, due to the correlation between the 
production and propagation mechanisms of secondary antideuterons and of the secondary 
nuclei on which the propagation models are tuned. However, for small kinetic energies, 
the propagation history of antideuterons, which is dominated by the adiabatic energy loss, 
is very different to that of other secondary nuclei. This contribution is proportional to 
Vc (see Eq. (10)), i.e. a larger value of the convection velocity implies a larger adiabatic 
energy loss. Therefore, the background flux is significantly lower for the MAX propagation 
parameters, as in this case the convection velocity Vc = 5km/s is rather low, while for the 
MED and MIN values this parameter is larger and nearly the same, with Vc = 12km/s 
{Vc = 13.5km/s) for MED (MIN). This source of uncertainty, nonetheless, is expected to be 
reduced in the near future with the improved measurements by AMS-02 of the flux ratios of 
stable secondary-to-primary and unstable cosmic ray species, such as boron-to-carbon and 
beryllium-lO-to-beryllium-9 |59j . 

A further source of uncertainty stems from the modeling of the non-annihilating rescat- 
tering cross sections, which enters the computation of the tertiary source term defined in 
Eq. (14). As explained above, there is scarce experimental information on the cross sections 
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Kiuetic energy per nucleon [GeV/u] Kinetic energy per nucleon [GeV/n] 

Figure 8: Impact of various sources of uncertainty on the calculation of the cosmic an- 
tideuteron flux: value of the coalescence momentum, po (upper left), choice of the propaga- 
tion parameters (upper right), size of the non-annihilating antideuteron-proton rescattering 
cross section (lower left) and modeling of the non-annihilating rescattering cross section 
(lower right). See the text for details. 



of the rescattering processes d + p/He d + X, therefore it is necessary to make several 
assumptions and estimations in their evaluation, making a definite error estimation very 
difficult. We investigate in Fig. [8| lower left plot, the impact on the result when doubling 
or halving our estimation for cr^p"'^""^, which amounts to a pretty conservative error of at 
most ~ 30%. This rather mild dependence is connected to the fact that the dominant 
effect governing the low-energy part of the antideuteron background flux is the adiabatic 
energy loss and not the non-annihilating rescattering. Lastly, we compare the result using 
the method described in this section with that followed in [601 H9] based on the "limiting 
fragmentation hypothesis" [51], and which assumes that the differential cross section for the 
NAR process is a simple step function with respect to the incident kinetic energy of the 
antideuteron [TSl [IS] : 

da^-g {dm) + p -> d{T,) +X) _ atot {d{TL) +p^d + X) 

d7> = ^^^^ ^^^^ 

The results are shown in Fig. [8} lower right plot, and are practically the same for both 
parametrizations. Again, this source of uncertainty could be better understood and reduced 
with dedicated experiments studying antideuteron-nucleus collisions. 
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4 Conclusions 

In this work we have investigated the production of cosmic antideuterons in high energy col- 
hsions of cosmic rays with the interstellar matter. An accurate calculation of the secondary 
flux is important not only to assess the prospects to observe antideuterons at AMS-02 or 
GAPS, but also to determine the background fluxes to search for primary antideuterons from 
exotic sources. This flux has been calculated in the past using a "modified factorized coales- 
cence model" which, due to the scarce experimental information on antideuteron production 
in high energy collisions, required a number of approximations and ad hoc assumptions. We 
have argued that some of the underlying assumptions of this approach are not validated 
by experimental data. Therefore, we have instead performed a full cvent-by-evcnt simu- 
lation of the antideuteron production implementing the coalescence model in the Monte 
Carlo generator DPMJET-III, which was designed to simulate low-energetic hadron-hadron 
collisions, such as those relevant for experiments searching for cosmic antideuterons. Then, 
we used an appropriate diffusion-convection equation for charged cosmic rays in order to 
translate the secondary antideuteron source spectrum into the expected flux of background 
antideuterons at Earth. We find a secondary antideuteron flux approximately a factor of 
two smaller than the most recent analyses and which amounts to 0.13 background events 
at AMS-02 and 0.024 events at GAPS ULDB. Lastly, we have also studied various sources 
of uncertainties affecting the computation of the antideuteron background flux and we have 
emphasized the importance of dedicated experiments analyzing antideuteron production in 
proton-proton collisions at y/s ~ 10 GeV to reduce the uncertainties. 
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